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We carry out high-precision calculation of parity violation in cesium atom, reducing theoretical 
uncertainty by a factor of two compared to previous evaluations. We combine previous measurements 
with our calculations and extract the weak charge of the '^^^Cs nucleus, Qw = — 73.16(29)exp(20)th- 
The result is in agreement with the Standard Model (SM) of elementary particles. This is the 
most accurate to-date test of the low-energy electroweak sector of the SM. In combination with the 
results of high-energy collider experiments, we confirm the energy- dependence (or "running") of the 
electroweak force over an energy range spanning four orders of magnitude (from ~ 10 MeV to ~ 100 
GeV). Additionally, our result places constraints on a variety of new physics scenarios beyond the 
SM. In particular, we increase the lower limit on the masses of extra Z-bosons predicted by models 
of grand unification and string theories. 
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Atomic parity violation places powerful constraints on 
new physics beyond the Standard Model (SM) of ele- 
mentary particles [l|, The measurements are inter- 
preted in terms of the nuclear weak charge Qw, quan- 
tifying the strength of the electroweak coupling between 
atomic electrons and quarks of the nucleus. Here we re- 
port the most accurate to-date determination of this cou- 
pling strength by combining previous measurements 0, 0| 
with our high-precision calculations in cesium atom. The 
result, Qvv("3Cs) = -73.16(29)cxp(20)th, is in a perfect 
agreement with the prediction of the SM. In combina- 
tion with the results of high-energy collider experiments, 
our work confirms the predicted energy-dependence (or 
"running" ) of the electroweak interaction over an energy 
range spanning four orders of magnitude (from ^ 10 
MeV to ~ 100 GeV). The attained precision is impor- 
tant for probing "new physics". As an illustration, we 
constrain the mass of the so-far elusive particle - the 
extra Z boson {Z'). Z' are hypothesized to be carriers 
of the "fifth force" of Nature, and they are abundant 
in models of grand unification and string theories 
In particular, SO(IO) unification predicts a Z'-boson de- 
noted as Z'-^. A direct search at Tevatron collider Q 
yielded Mz^ > 0.82TeV/c^. Our precision result implies 
a more stringent bound, Mz' > 1.3TeV/c^ If Z' is dis- 
covered at the Large Hadron Collider (LHC), where the 
mass scale reach is somewhat higher, our result would 
help in exacting Z' properties. 

Historically, atomic parity violation helped in estab- 
lishing the validity of the SM 0, i, H- While a num- 
ber of experiments have been carried out, the most ac- 
curate measurement is due to Wieman and collabora- 
tors \3\ . They determined a ratio of the parity noncon- 
serving (PNC) amplitude, -EpNC; to the vector transition 
polarizability, /3, Ep^c/P = 1.5935(56) mV/cm, on the 
parity-forbidden 681/2 7Sif2 transition in atomic Cs. 

The measurement, however, does not directly translate 



into an electroweak observable of the same accuracy, as 
the interpretation of the experiment requires input from 
atomic theory. The theory links Qw to the measured 
signal. In computations, Qw is treated as a parameter, 
and by combining computed -EpNC with measurements, 
the value of Qw is derived. The inferred Qw com- 
pared with the predicted SM value, either revealing or 
constraining new physics. So far the atomic-theory error 
has been a limiting factor in this interpretation. Here we 
report reducing this error, leading to an improved test of 
the SM. 
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Here D and ifw are electric-dipole and weak interac- 
tion operators, and Ei are atomic energy levels. In the 
electronic sector, the effective weak interaction averaged 
over quarks reads Hw = — ^QvK75p(r), where G_f is 
the Fermi constant, 75 is the Dirac matrix, and p{r) is the 
neutron-density distribution. ^^"^Cs nucleus has Z = 55 
protons and N — 78 neutrons. The value of Qw is given 
approximately by —N. 

Interpretation of the PNC measurements requires eval- 
uating Eq.([T]). Although the underlying theory of quan- 
tum electrodynamics (QED) is well established, the 
atomic many-body problem is intractable. Reaching the- 
oretical accuracy equal to or better than the experimen- 
tal accuracy of 0.35% has been a challenging task (see 
Fig. [T]). An important 1% accuracy milestone has been 
reached by the Novosibirsk [Ifl] and Notre Dame [llj 
groups in the late 1980s. More recently, several groups 
have contributed to understanding sub-1% corrections, 
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primarily due to the Breit ( mag netic) interaction and ra- 

[3 [ll [111 (reviewed in 



12, U 



diative QED processes 

[13] )■ The results of these calculations are summarized 
by the "World average '05" point of Fig. [1] which has a 
0.5% error bar reflecting this progress. As of 2005, the 
sensitivity to new physics has been limited by the accu- 
racy of solving the basic correlation problem. Here we 
report an important progress in solving it. 
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FIG. 1: (Color online) Progress in evaluating the PNC ampli- 
tude. Points marked Paris '86, Novosibirsk '89, Notre Dame 
'90 correspond to Refs. [Tl], [l^l, and Point "World av- 
erage '05" is due to efforts of several groups [l^. [l^, [T3. fisl, [T^ 
on sub-1% Breit, QED, and neutron-skin corrections reviewed 
in Ref. The strip corresponds to a combination of the 

Standard Model Qw with measurements [3, The edges of 
the strip correspond to zba of the measurement. Here we ex- 
press -BpNC in conventional units of i|e|aB {—Qw/N) x 10~^^, 
where e is the elementary charge and ob is the Bohr radius. 
These units factor out a ratio of Qw to its approximate value, 
-iV. 



We wish to evaluate accurately the sum ([T|). To this 
end, we solve the Schrodinger equation H\'^y) — Eyl'i'y) 
and find atomic wave functions and energies. Even in 
classical mechanics, the simpler three-hody problem can- 
not be solved in closed form. For Cs atom, one solves for a 
correlated motion of 55 electrons. The problem is simpli- 
fied by the fact that this atom has one loosely-bound va- 
lence electron v outside a stiff closed-shell core. Because 
of that, the problem can be efficiently treated within the 
many-body perturbation theory p^ . In this treatment, 
the exact many-body state \'^v) which stems from the 
approximate (Dirac-Fock) state is parameterized 

as I 'ft,) = Q l^'i"'') , where the many-body operator fl is 
yet to be found. It is expanded into a hierarchy of single-, 
double-, triple-, and higher-rank n-fold excitations. For 
example, double excitations (or simply doubles) result 
from a simultaneous scattering of two core electrons by 
their mutual Coulomb repulsion. 

Notice that for the 55 electrons of Cs this treatment 
would be exact by including 55-fold excitations. How- 
ever, manipulating such wave functions is impractical: 
for a basis set of 100 orbitals, one would require more 
than 100^^ memory units. This number exceeds the es- 
timated number of atoms in the Universe. Fortunately, 
contributions of high-rank excitations are strongly sup- 
pressed. Previous many-body calculations [Til. [20 | in Cs 



stored only single and double excitations. Already at this 
level the attained accuracy for the atomic properties of 
Cs was at the level of 1% or better. To systematically 
improve the accuracy, here we take advantage of modern 
computing resources and additionally store and manip- 
ulate triple excitations. This is a substantial step. For 
example, previous calculations [Tll| used less than 100 Mb 
of storage, whereas our calculations required 100 Gb; this 
is a factor of 1,000 increase in computational complexity. 

Our specific scheme [13, Hi], IH, [2^ of solving the 
atomic many-body problem is rooted in the coupled- 
cluster method [l9|. We refer to our approximation 
as the CCSDvT scheme (Coupled-Cluster approximation 
including Singles, Doubles, and valence Triples). Details 
will be provided elsewhere. The solution is ab initio rela- 
tivistic, as near the Cs nucleus (where the weak interac- 
tion occurs) the electrons move with speeds approaching 
the speed of light. To minimize human errors, the CCS- 
DvT code was developed independently by at least two 
persons. Complex derivations and coding were aided by 
symbolic algebra tools. An important proof of the code 
was made by computing properties of lithium atom [23[ . 
This atom has 3 electrons, making the CCSDvT approx- 
imation exact. We found in Ref. [1^ that experimental 
data for Li were reproduced numerically with an accuracy 
reaching 0.01%. 

Now we proceed to evaluating the PNC amplitude, 
Eq. II]), by directly summing over the intermediate nPi/2 
states {ll| . This implies computing wave functions and 
energies of the 6S'i/2, 7S'i/2, and nPi/2 states, forming 
matrix elements, and substituting them into Eq. We 
employ a computationally expensive CCSDvT method 
only for matrix elements involving rt = 6, 7, 8, 9 ("main" 
term) and compute suppressed contributions of rt > 10 
and core-excited states ("tail" term) with less accurate 
methods. 

Our results for the PNC amplitude are presented in 
Table IT] The upper panel of the table lists contribu- 
tions due to the Coulomb interaction of electrons with 
the nucleus and other electrons. The lower panel sum- 
marizes well-established non-Coulomb contributions such 
as Breit, radiative (QED), and other smaller corrections. 
Estimated uncertainties are listed in parentheses. 

We start by assessing the accuracy of the employed 
CCSDvT approximation. This determines uncertainty 
of the "main" term contributing 99% of SpNC- Prop- 
erties of low-energy states have previously been mea- 
sured, and we quantify theory uncertainties by compar- 
ing these data with our ab initio results. For consistency 
we add QED, Breit, and nuclear-structure corrections to 
our Coulomb-correlated results. We find that the exper- 
imental energies are reproduced with an accuracy of 0.1- 
0.3%. Dipole matrix elements enter the PNC amplitude 
directly and are derived from atomic lifetime measure- 
ments. Relevant dipoles are compared in the lower panel 
of Fig.O the CCSDvT values are within the error bars of 
the experiments. Finally, since the hyperfine constants A 
arise due to interactions of electrons with nuclear mag- 
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TABLE I: Contributions to the parity violating amplitude 
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75'i/2 transition in Cs in units of 



Coulomb interaction 




Main (n = 6 - 9) 


0.8823(18) 


Tail 


0.0175(18) 


Total correlated 


0.8998(25) 


Corrections 




Breit, Ref. \12} 


-0.0054(5) 


QED, Ref. \16} 


-0.0024(3) 


Neutron skin, Ref. [1^ 


-0.0017(5) 


e — e weak interaction, Ref. [11] 


0.0003 


Final 


0.8906(26) 



smaller non-Coulomb corrections and its uncertainty was 
estimated by adding individual uncertainties in quadra- 
ture. Previous calculations [ll|j report values larger 
by 0.9% than our 0.27%-accurate result. The difference 
is due to our inclusion of additional many-body effects, 
shown in Fig. [31 Direct contribution of triple excitations 
to matrix elements accounts for a 0.3% shift and dress- 
ing of matrix elements for another 0.3%. The remaining 
0.3% comes from a consistent removal of QED and Breit 
corrections from experimental energies during the semi- 
empirical fit. 



netic moments, matrix elements of the weak interaction 
{nSi/2\Hw\n' P1/2) may be tested by forming the geo- 
metric mean ^^nSi/a^"' P1/2 > l2fl]- Deviations of 
these combinations from experimental data are shown in 
the upper panel of Fig. [2l We find that the standard 
deviation of theoretical values from experiment is 0.2%. 
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FIG. 2: (Color online) Deviations of computed values (red 
filled circles) from experimental data (centered at zero). The 
upper panel displays combinations of magnetic hyperfine 

which mimic matrix el- 



structure constants ^ Ans^^/ 

ements of the weak interaction. For these combinations, ex- 
perimental error bars are negligible compared to the theo- 
retical accuracy. The lower panel exhibits deviations of the 
computed dipole matrix elements from the most accurate ex- 
perimental results [2^. l25j. 



Overall agreement of theoretical data with experiments 
indicates that the average accuracy of the CCSDvT ap- 
proximation is 0.2% and we assign an error of 0.2% to 
the main term. Additionally, our semi-empirical fitting 
to experimental energies modifies the main term by 0.2%, 
which is consistent with the above error estimate. Fi- 
nally, the "tail" was computed using a blend of many- 
body approximations and we assign a 10% uncertainty to 
this contribution based on the spread of its value in dif- 
ferent approximations. The final result (Table |T| includes 




FIG. 3: Many-body diagrams responsible for the shift of 
the PNC amplitude compared to previous calculations. Top 
row; sample direct contributions of valence triples to ma- 
trix elements (wavy capped line) [22|. Bottom row: iterative 
equation for line dressing of the hole line in expressions for 
matrix elements [U (similar equation holds for particle lines; 
exchange diagrams are not shown). 



With the computed E'pnc we proceed to extracting 
the electroweak observable. The experiment [3j deter- 
mined the ratio Ep^c/P = 1-5935(56) mV/cm, /? being 
the vector transition polarizability. The most accurate 
value of P comes from a combined determination 0, [20| , 
(3 = -26.957(51)a|. With this /3, we arrive at the nu- 
clear weak charge 



^H-C^'Cs) = -73.16(29)exp(20)th 



(2) 



where the first uncertainty is experimental and the sec- 
ond uncertainty is theoretical. Taking a weighted av- 
erage, (3 = — 26.99(5)a^, of two determinations [1, [25j 
results in Qwi^^^Cs) = -73.25(29)cxp(20)th. Both val- 
ues are in a perfect agreement with the prediction of the 
SM, = -73.16(3) of Ref. [H. 

Our result plays a unique, and at the same time com- 
plementary, role to collider experiments. For ^^'^Cs atom 
the relevant momentum transfer is just ~ 30 MeV [lH . 
but the exquisite accuracy of the interpretation probes 
minute contributions of the sea of virtual (including so- 
far undiscovered) particles at a much higher mass scale. 
The new physics brought by the virtual sea is phe- 
nomenologically described by weak isospin-conserving S 
and isospin-breaking T parameters j27j : AQw — Qw — 
qSM = -0.800 - 0.007T. At the la-level, our result 
implies \S\ < 0.45. Parameter S is important, for ex- 
ample, in indirectly constraining the mass of the Higgs 
particle [131 ■ Similarly, the extra Z boson, Z'^, dis- 
cussed in the introduction, would lead to a deviation [l[ 
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AQw w 84{Mw/Mz^)'^, where Mw is the mass of the 
W boson. Our result implies 50% chance that there is 
Z' (i.e., AQw > 0). We find (at 84% confidence level, 
including M^/^ = oo) Mz^ > 1.3TeV/c^, raising the 
present lower bound on the mass from direct collider 
searches @ . Our raised bound on the Z' mass carves out 
a lower-energy part of the discovery reach of the LHC. 

Our result confirms fundamental "running" (energy- 
dependence) of the electroweak force [H, [23] ■ The in- 
teraction strength of particles depends on their relative 
collision energy E: at higher energies the collision part- 
ners tend to penetrate deeper inside the shielding clouds 
of virtual particles surrounding the particles. According 
to the SM, the interaction strength at low energies dif- 
fers by about 3% from its measured value at 100 GeV. 
Compared to collider experiments, our result provides a 
reference point for the least energetic electroweak interac- 
tions. Notice that the previous analyses 17] of PNC 
in Cs were consistent with no running [30|. With our 
Qw, we find the effective interaction strength (we use 
scheme of Ref. jsl), sm'^9'i^{E ^ 0) = 0.2381(11), 0w 
being the Weinberg angle. The uncertainty is somewhat 



better than that of the previous most precise low-energy 
test of the electroweak sector obtained in the electron 
scattering experiment at SLAC [S^l- Our result is in 
agreement with the SM value [sT] of 0.2381(6). While an 
earlier evidence for running of siia^Ow has been obtained 
at SLAC [121, the prediction of the SM was outside of 
their error bars. In this regard, in addition to placing 
important constraints on new physics beyond the SM, 
our work provides a higher-confidence confirmation of the 
predicted running of the electroweak coupling at low en- 
ergies. In combination with the results of high-energy 
experiments at SLAC and CERN [2^ , our work confirms 
the predicted running of the electroweak interaction over 
an energy range spanning four orders of magnitude (from 
- 10 MeV to - 100 GeV). 
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